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^ . Abstract 

Complete topological classification of solutions in 50(3) symmetric 

^ \ Ginzburg-Landau free energy has been performed and a new class of solutions 

in weak external magnetic field carrying two units of magnetic flux has been 

^ identified. These solutions, magnetic skyrmions, do not have singular core 

Tlj- ' like Abrikosov vortices and at low magnetic field become lighter for strongly 
(N ■ 

^ ^ type n superconductors. As a consequence, the lower critical magnetic field 

I Hci is reduced by a factor of logK. Magnetic skyrmions repel each other as 

I 1/r at distances much larger then magnetic penetration depth A forming rela- 

^ ^ ^ tively robust triangular lattice. Magnetization near Hd increases gradually as 

Q I (H — Hci)^. This behavior agrees very well with experiments on heavy fermion 

superconductor UPt^. Newly discovered Ru based compounds Sr2RuOi and 

^ ' Sr2Y Rui^xCuxOq are other possible candidates to possess skyrmion lattices. 

Deviations from exact 50(3) symmetry are also studied. 
PACS numbers: 74.20.De, 74.25.Ha, 74.60.Ec, 74.70.Tx 
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I. INTRODUCTION 



Rich variety of novel magnetic properties can be found in superconductors with uncon- 
ventional type of pairing symmetry. At present several examples of unconventional super- 
conductors are known. The first of them is, of course, a family of Tc cuprates. In connection 
with them a lot of recent efforts were devoted to study of d - wave pairing of various types 
with possible admixture of s - wave. On the other hand, a triplet type of pairing is believed 
to exist in UPt^ and in some other heavy fermion compounds [Q. It is also suspected 
to occur in recently discovered new classes of Ru based superconductors, layered perovskite 
Sr2RuOi @,|] and bulk compound Sr2Y Rui-xCuxO^ ^ which has a double perovskite 
structure. Since all mentioned superconductors are of strongly type II, vortices play the 
major role in their thermodynamical properties. In high superconductors, despite fun- 
damental differences in mechanism and microscopic properties compared to conventional 
superconductors, vortices are quite similar to conventional Abrikosov vortices. The reason 
is that there exists a dominant single order parameter field: d - wave condensate. Small 
(sometimes quite important) deviations can be accounted for due to admixture of the s - 
wave component. Then, the order parameter is effectively multicomponent. This property 
leads generally to various new effects like nonaxisymmetric vortices and phase tran- 
sitions within flux line lattices near H(.2 §]• Similar phenomena exist and should be even 
more pronounced in the systems with intrinsically multicomponent superconducting order 
parameter |1T0| , [TT|| like heavy fermions compounds. 

The situation in triplet superconductors might be more exotic. The order parameter is 
necessarily multicomponent. In addition, under certain conditions the rotational symmetry 
(at least approximate one) between different components might exist. In that case vortices 
are not the only type of topological solitons which can carry magnetic flux through the sam- 
ple. The corresponding phenomenological Ginzburg-Landau theory has the order parameter 
of a vector type with a continuous symmetry. It is known in the theory of superfiuid ^He 



that in such a system there exist topological defects which have no singularities even 
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within the London approximation. On the other hand, vortices have a singularity at their 
core, at least for one of components of the order parameter. This makes their energy roughly 
proportional to log/t, similarly to the case of a standard Abrikosov vortex. Therefore, for 
sufficiently large n vortices are expected to be heavier then nonsingular topological defects 
and the latter become most likely candidates for thermodynamically stable configurations 
of the order parameter field into which a homogeneous superconducting (Meissner) state 
transforms under the action of an external magnetic field. 

It is the purpose of the present paper to investigate this possibility in detail. We find 
such a nonsingular in the London approximation solution, the magnetic skyrmion, describe 
its structure and show that it is energetically favorable over Abrikosov vortex in wide range 
of Ginzburg-Landau parameter n values. Lattices of magnetic skyrmions are particularly 
important at fields near the lower critical field. Most striking effects are reduction of Hd by 
a factor of logn and dramatic change in the behavior of magnetization near Hd. We also 
investigate what happens to magnetic skyrmions when the continuous symmetry breaking 
terms are introduced into the free energy. It is shown they survive under small perturbations 
and gradually evolve to another still nonsingular configurations under large perturbations 
of a certain type. 

Magnetic skyrmion lattices may have been already experimentally observed in UPt^. 
Magnetization curves of near Hd [|l3l , p!4| are rather unusual (see Fig.l of |15| in which 



a short account of this work was presented). Theoretically, if the magnetization is due 
to penetration of vortices into a superconducting sample then one expects —AttM to drop 
with an infinite derivative at Hd- On the other hand experimentally —4ttM continues to 
increase smoothly. Such a behavior was attributed to strong fiux pinning or surface effects 
13| . However both experimental curves in Fig.l of [|T5[, as well as the other ones found in 



literature, are close to each other if plotted in units of Hd- We propose a more fundamental 
explanation of the universal smooth magnetization curve near Hd- If one assumes that 
fluxons are of unconventional type for which interaction is long range then precisely this 
type of magnetization curve is obtained. Indeed magnetization near Hd due to fluxons 



carrying units of flux $o = hc/2e, with line energy e and mutual interaction V{r), is 
found by minimizing the Gibbs energy of a very sparse triangular lattice: 

where oa = ($o/-Bv^)^ is lattice spacing. When V{r) ~ exp[— Ar], the magnetic induction 
has the conventional behavior B ~ [log {H — Hci)]^"^ ||16[, while if it is long range, V{r) ~ 
1/r", then one finds B ~ [H — Hd)^^^ ■ The physical reason for this different behavior is 
very clear. For a short range repulsion, if one fiuxon penetrated the sample, many more 
can penetrate almost with no additional cost of energy. This leads to the infinite derivative 
of magnetization. On the other hand for a long range interaction making a place for each 
additional fiuxon becomes energy consuming. Derivative of magnetization thus becomes 
finite. 

The remainder of the paper is organized as follows. In section II we present the 5*0(3) 
symmetric model and note that it is an excellent approximation to certain successful models 
of UPts |T^ as well as to others ||T^. The London approximation is developed. In section 
III we perform complete topological classification of solutions and find that the magnetic 
skyrmion carries two units of magnetic fiux. General form of cylindrically symmetrical 
solutions is given. In section IV we determine magnetic skyrmion lattice structure, Hd 
and the magnetization curve. An example of deviations from exact 5*0(3) symmetry is 
considered in section V. More specifically, we address the case of Zeeman like interaction 
relevant to Sr2Y Rui^^CuxO^ [|1^ system which initially motivated us to search for exotic 
vortices. Section VI contains discussion of the results and possibilities to experimentally 
observe various effects of magnetic skyrmions. 



II. THE MODEL 
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A. Ginzburg-Landau free energy functional and its symmetries 



Let us consider a model Ginzburg-Landau theory with the order parameter ipi{r) being 
a three dimensional {i = 1, 2, 3) complex vector. It is convenient to consider index i as a 
spin in the case of weak spin-orbit coupling in the pairing channel, but this is not necessary 
interpretation. The case of strong spin - orbit interaction can also be addressed provided 
some modifications of the free energy functional are made. Average spin of the Cooper pair 
at a specific point in the space is given by 

S,{f) = rj{^Heijk)Mr)- (2) 

The material under study is assumed to be isotropic. Extensions of our results to anisotropic 
situations are discussed in the section V. 

The Ginzburg - Landau free energy functional of the system has the form: 

F,ot = -a^^^: + ^(^.^n' + y i^^^.r, (4) 

F,ra, = £-^{V,ij,)(V,i:,r, (5) 

where = dj—i{e*/hc)Aj are covariant derivatives, Bj = {VxA)j, m* > is effective mass 
of the pair and e* is effective charge of the pair. For a superconducting phase to exist the 
coefficient a should be positive below the phase transition point and we set a = a'{Tc — T) 
with a' > 0, while for positive definiteness of the potential the other coefficients of eq.(^ 
should satisfy /?i > and P2 > —Pi- 

The free energy density eq.@ has the following independent symmetries. The spin 
rotations, forming a group S'0spm(3), act on the index i of the order parameter field, so that 
it transforms as a vector. Two dimensional (orbital) space rotations, forming a different 
SOorbiti'^) group, act on spatial coordinates Xj and the electric charge transformations, 
forming a U{1) group, rotate the complex phase of the order parameter. Note that in eq.(^ 
we assumed that external magnetic field H is oriented along z direction. We will consider 
only configurations invariant under translations in that direction or the thin film geometry. 
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First we consider the case of zero external magnetic field. Pure superconducting (Meiss- 
ner) phases that appear below are found by minimization of Fpot with respect to ip*. This 
is conveniently done making use of the following parametrization of the order parameter 
vector: 

'j/' = V'jej = /(ncos + i m sin 0), (6) 

where />O,O<0<7r/2,n and m are unit vectors that are arbitrarily oriented with 
respect to some fixed coordinate system in the spin space with orthonormal basis 61,62,63. 
There exist two phases, depending on the sign of the coefficient (32 '■ 

I:/?2>0 V^=/^±p, n±m,0 = 7r/4, f = j^. (7) 

II:/32<0 i^=fe'^n, n = ±m- f = (8) 

Pi + P2 

Combining eq.(|D and eq.(^ one obtains S = p sin 20 Z, where I = n x m. In the phase I 
the projection of the spin of a Cooper pair S on the vector I is equal to either +1 or — 1, 
refiecting spontaneous time reversal symmetry breaking. In the phase II this projection is 
always zero. 

B. Theories of triplet superconductors and terms breaking S'Ospin(3) symmetry 

Obviously, the model of the previous subsection is an idealization of the actual situation 
in triplet superconductors. In this subsection we note that some successful models of f/Pta, 



notably that of Machida et al differ from this model only by small less symmetric 
terms. These terms could be considered as small perturbations, at least in some regions 
of if — T diagram. The asjTiimetries are of several types. First, the space symmetry 
SOorbiti'^) is normally broken down to some crystallographic point group of a given material 
{Dgh for UPts, D^h for Sr2Ru04, D2h for Sr2Y Rui^xCuxO^) . The effective mass m* then 
becomes a symmetric tensor m*^. Second, the spin S can be coupled to the magnetic field. 
This explicitly breaks SOspin{,i) down to SOspin{'^)- Separate spin and orbital symmetries 
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S0spin{'2) ® SOoThitiX) are broken down to diagonal SOtotC^) as well. The various types of 
perturbations are 

AFpot2 = a' [(T, - T^''Mi\' + (T, - , (9a) 

AFpotA = PsilH' - l^iT - 1^21')' + - l^iT - 1^21'), (9b) 

AFgrad = K [(P,^,)(PiVj)* + {ViiJi){V,ijjy] , (9c) 

AFzeeman = fiS ■ B, (9d) 

AFnonlin = Ax\^f ' (9e) 

We estimate their coefficients for the case of UPt^. Some models of UPt^ do not have 
three dimensional complex order parameter and therefore will not be addressed here. Exam- 
ples include Eig singlet pairing [^], i?2gtriplet pairing accidental degenerate AB model 



. On the other hand, the model of weak spin-orbit coupling developed by Machida et al 
p!7| is of the type we are interested in. In this model asymmetric terms are very important 
in explaining the double superconducting phase transition at zero external magnetic field. 
However, they are small in the low temperature superconducting phase (phase B) well below 
its critical temperature T ^ ~ AbK and at low magnetic fields H ~ Hd- Indeed, for 
the quadratic terms one gets from experiment '^"7^"^ ~ .2, ^^i'^^^ ^ _q5_ ^jj-^g quartic terms 
eq.(^) are order of magnitude smaller. The corrections to the gradient terms are very small 
K/ (2^7^) ~ -^^ ^^"^ /^/ (2m*) ~ '^^^ coefficient of the nonlinear coupling term eq.(|9cD 



is negligible: {^^Hl^ / — 10"^. Another model of UPt^, which has similar structure 
to eq.(|]), is the accidental degenerate AE |T^. Estimates are similar with exception of 
K/ (2^) "which is now of order one (found from ffiting transition lines near Hc2)- 

In general, topological solitons exist even in those cases when the symmetry is weakly 
broken. In section V we consider in detail influence of one symmetry breaking term, Zeeman 
coupling, eq. (pBj), in connection to a new material Sr2Y Rui^xCuxOq. We show that it do 
not affect stability of solitons that we investigate in this paper 
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C. London approximation 



The London approximation assumes that the order parameter has the form determined 
by the potential part of the free energy eq.(^. In particular, the modulus of the order 
parameter is fixed. Any variations of the order parameter field over the space are only due 
to changes of the degeneracy parameters which parametrize the vacuum manifold. From 
this standpoint in usual s-wave superconductors there are no topological solitons within 
the London approximation. Famous Abrikosov vortex has a core -a region in which the 
modulus of the order parameter varies significantly and vanishes at some point. A vortex 
can be incorporated into the London approximation at the cost of singularities: vortex core 
is assumed to shrink to a point in which energy diverges logarithmically. Accordingly, a 
cutoff, the correlation length, should be introduced and one obtains log k dependence for a 
vortex line tension. As discussed above, this means that if there exists a nonsingular solution 
it is bound to become energetically favorable for k large enough. 

Below we concentrate on the properties of nonunitary phase I, eq.(0), of a triplet super- 
conductor near the lower critical field Hd- This phase is always assumed when we refer to 



the superconducting state. We define magnetic penetration depth A = j^\y^4^, coherence 



length ^ = 7i j \J lam* ^ flux quantum $o = hcje* and Ginzburg-Landau parameter k = X/S,- 
For convenience of the following discussion we express all physical quantities in dimensionless 
units as follows: 



..^.,F^£,P, P.fr-, A.^^a, B.^, (10) 



The "tilde" marks will be omitted hereafter. 

In order to determine the degeneracy parameters we consider the symmetry breaking 
pattern of the superconducting state. Both the spin rotation SOgpini^) symmetry and the 
superconducting phase U{1) symmetry are spontaneously broken, but a diagonal subgroup 
f/(l) survives. It consists of combined transformations: rotations by angle i) around the axis 
/ which are accompanied by gauge transformations e*'^. These combined transformations 
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together with rotations of vector / itself, form the vacuum manifold. The vacuum manifold 
is isomorphic to the 50(3) group. Our aim is to find nonsingular topological line defects in 
this case. We choose a triad of orthonormal vectors n, m, I to be the degeneracy parameters. 
From the definition of these vectors the following important relations can be derived: 

ndiTn = —mdifi, (11) 
{d,nf + (d.my = 2{nd,mf + (9,/?, (12) 
£pqslp{dilq){djls) = {dinp){djmp) - {dimp){djnp). (13) 

To obtain the free energy density of the London approximation we substitute ip in the 
form eq.(|^ into the gradient part, eq.(|^), of the total free energy functional and make use 



of eg. ([Tl|) and eq.(|T^). After some algebra we get: 

Fl = ^ {d,f)' + {nd,m - aif + bl (14) 

Varying energy functional with respect to vector potential a one obtains the supercurrent 
equation: 



UpViTLp — a = V X (^V X aj = j, (15) 

where the Maxwell equation was used. Eq.(|T5[) shows that the superconducting velocity (in 
units of h/m*) is given by 

UpVmp = -V^. (16) 

Thus, the angle i), which specifies the position the pair of perpendicular unit vectors n and 
m in the plane normal to vector /, takes the role of superconducting phase in the present 
case (see Fig.l). Other field equations are most easily obtained by considering Fl{1, n, m) as 
a functional of I and n only and performing conditional variation with constraints /■ = 0, 
P = = 1. This procedure yields independent equation for /: 

Af- Al) + 2jkilx dkl) = 0. (17) 
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III. TOPOLOGICAL CLASSIFICATION OF SOLUTIONS IN LONDON 



APPROXIMATION 

In this subsection we develop a classification scheme for the finite energy solutions to 
our model in the London approximation derived above. The main results is that the Lon- 
don equation, eq.(p!5D and eq. ([T7|) , in the presence of the magnetic flux admit nonsingular 
topologically stable solutions. This class of solutions contains cylindrically symmetric ones. 



A. General topological analysis 



Let us consider boundary conditions for a superconductor which extends over the whole 
space. The free energy density eq.(0) is positive definite and contains term. It follows 
that magnetic field vanishes at spatial infinity. Then one has to specify the triad n, fh, I at 
different distant points. The corresponding (first) homotopy group of vacuum manifold is 



7Ti{SO{3)) = Z2 |12|. It yields a classification of finite energy solutions into two topologi- 
cally distinct classes. This classification is too weak, however, because it doesn't guarantee 
nontrivial flux penetrating the plane. We will see that conflgurations having both "parities" 
are of interest. 

In the presence of the magnetic flux, the conflgurations are further constrained due to 
the flux quantization condition. The vacuum manifold is naturally divided into 5*0(3) — *• 
5*0(2) 5*2, where the 5*2 is the direction of / and the 5*0(2) is the superconducting phase 
deflned in eq.(|l^). For given number of flux quanta N = $/$05 the phase makes 
N winds at inflnity, see Fig.2. The flrst homotopy group of this part is therefore flxed: 
TTi (5*0(2)) = Z. If, in addition, vector / is flxed throughout the volume of a superconductor 
there is no way to avoid singularity in the phase ^9. It becomes ill deflned at some point 
and, accordingly, the modulus of the order parameter have to vanish there. Destruction of 
the superconducting state takes place in rather small area, especially for large k. Thus, we 
arrive at usual picture of Abrikosov vortex. 
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However, the general requirement that a solution has finite energy is much weaker. It 
tells us that the direction of I should be fixed only at infinity This follows from the presence 
of {d.d) term in Fl (see eq.(|14|)) which cannot be "gauged away" as the corresponding 
term for the 5*0(2) part. A relevant homotopy group is 7r2(S'2) = Z. The second homotopy 
group appears because constancy of / at infinity (say, up) effectively " compactifies" the 
two dimensional physical space into 5*2. One can have topologically nontrivial configuration, 
skyrmions, which are markedly different from vortices. Unit vector / can nontrivially wind 
towards the center of the texture. New topological number Q should be introduced [p^] : 



Q = ^ / ^iJ ^ {^i^^ djl) dS. (18) 

Configurations of the order parameter field with topological number Q = 1 have vector / 
fiipping its direction from up to down (or from down to up) until it reaches the center of 
the texture from an infinitely remote point (see Fig.2). 

To summarize, configurations fall into classes characterized by two integers and Q. 
The "parity" of the more general topological analysis is just Q = A^(mod 2). Due to the pres- 
ence of two topological numbers an interesting possibility arises. There exists topologically 
nontrivial configuration that preserves the modulus of the order parameter (see eq.(|^)) at 
every point. We call these regular solutions magnetic skyrmions. For them these two topo- 
logical numbers are related to each other. We find this relation integrating the supercurrent 
equation, eq.(|15D along a remote contour and using of the identity eq. ([13|) : 

Q = N/2. (19) 

The lowest energy solution within the London approximation corresponds to A^/2 = Q = 

±1. 



B. Cylindrically symmetric magnetic skyrmions. 



In the class of solution A^/2 = Q = —1 there are ones possessing cylindrical symmetry. 
We will describe them in polar coordinates, p and (p. The triad n, m, / has form: 
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I = Cz COS 0(p) + Cp sin 9(p), 

n = (6*2 sin 0(p) — Cp cos 0(p)) sin + cos ip (20) 
m = (e^ sin 0(p) — Cp cos 0(p)) cos ip — sin y), 

where G is the azimuthal angle of / (see Fig.l). This choice corresponds to the situation 
when the pair of perpendicular vectors n and m winds twice as a distant circle on Fig. 2 is 
completed. Due to cylindrical symmetry of the solution in question function 0(p) satisfies 
boundary conditions = vr at p = and B = at p ^ oo. 

The free energy of the magnetic skyrmion per unit length takes form: 

Cms ~l~ ^cur ~l~ ^mag (^-^) 



£s = J pdp 



1 dQ\ sin^e 

2 [l^J ^ 2p2 

\ 2 



(22) 



, / 1 + cos , 
pdp + a (23) 



P 

?2 



pdpB^ =/pdpg + g . (24) 

where energy is measured in units of eg = • The first part e^s is the same as in standard 

nonlinear a— model without magnetic field [^. The second term Scur is analogous to the 
supercurrent contribution in the London approximation of the usual superconductor [jl6 



The third term is the magnetic energy. Eq. (|2l|) shows that a singularity at p = is absent 
(integrand converges) since 1 + cos 0(0) = 0. 

Actual distribution of magnetic field and order parameter in this case can be found from 
the following system of equations: 

0" + le' = -^f^±^ + 2aV (25) 
P P \ P J 

a" + --^-a = -(1 + COS0). (26) 
p p2 p 

In the next section we solve this equation. 
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IV. MAGNETIC SKYRMION SOLUTION 



A. Blow up of single skyrmion by magnetic field 



The general form of the solution of eq.(p5D and eg . (|2q) is given on Fig. 2. The orientation 
of the unit vector I (solid arrows) forms a skyrmion of 5*0(3) invariant a- model p3[. The 



phase t? makes two rounds at infinity (clock inside small circles on the "infinitely remote" 
circle). If magnetic field were absent there are infinitely many degenerate solutions 

e,(p) = 2arctan(5/p) (27) 

which have the same energy e = 2 for any size of the skyrmion S. The skyrmion of nonlinear 
a- model possesses a scale invariance. This degeneracy in various physical problem is lifted 
by perturbations. In some physical situations the skyrmion is stabilized by four derivative 
terms p3|, sometimes it shrinks and sometimes blows up. In the present context magnetic 
field lifts the degeneracy and we prove below that the skyrmion blows up. Of course is there 
are many skyrmions present their repulsion with stabilize the system. This is discussed in 
the next subsection. 

To prove that the skyrmion blows up, we explicitly construct variational configurations 
and show that as size of these configurations increases, the energy is reduced to a value 
arbitrarily close to the absolute minimum of Ems = 2. 

The first term in the energy eq.(PTD Es is the usual expression for the energy of the 
skyrmion. It is bound from below by the energy of usual skyrmion e = 2. To construct 
a variational configuration for 0, we pick up one of these solutions eq.(27) of certain size 
6. The second term Ecur, the "super current" contribution is positive definite. Therefore 
its minimum cannot be lower then zero. One still can maintain the zero value of this 
term when the field G is a skyrmion. Assuming this one gets the relation between a and 
6 : a{p) = — i+^i© = —-^£-^_ The magnetic field contribution (which is also positive 
definite) for such a vector potential is: Emag = afa • To sum up, the energy of the configuration 
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is e = 2 + gp-- It is clear that when 5 — > oo, we obtain energy arbitrarily close to the lower 
bound of e = 2. The skyrmion therefore blows up. 

We also solved eqs. (p5D-(p^) numerically on the segment of p from to a cutoff p^ax 
with boundary conditions h\p=p^^^ = and (9' + 0/p) \p=p^^^ = 0. The second boundary 
condition allows us to approach the correct asymptotic behavior of at infinity ~ 1/p 
which follows from eqs.(p5D-(p6|). The results for the distribution of magnetic field for Pmax 
ranging from 50 to 600 are presented in Fig. 3. One clearly sees that as the cutoff increases 
the magnetic field at the center p = decreases and the flux spreads out over larger area. 
This is in accord with the variational proof above. 



B. Skyrmion lattice and Hd 



Skyrmions repel each other, as we will see shortly, and therefore form a lattice. Since they 
are axially symmetric objects, the interaction is axially symmetric and hexagonal lattice is 
expected (see Fig. 4). Assume that lattice spacing is oa. At the boundaries of the hexagonal 
unit cells the angle is zero, while at the centers it is vr. Magnetic field b is continuous 
on the boundaries. Therefore, to analyze magnetic skyrmion lattice we should solve the 



equations (P5|)- (|26|) on the unit cell with such a boundary conditions demanding that two 



units of flux pass through the cell (by adjusting the value of magnetic field on the boundary). 

ol/4 

We approximate the hexagonal unit cell by a circle of radius R = ^7^=0 a having the same 



/2lT 

area. Fig. 4. 

We performed such calculations for R from R = 5 till R = 600 using finite elements 
method. The result is presented in Fig. 5. The energy per unit cell is described well in a 
wide range of R (deviation at i? = 10 is 1%) by an approximate expression 

eceii = 2 + ^ (28) 

The dominant constant contribution to the energy at large R comes like in the analytical 
variational state above from the first term Eg in the integrand of eq. (pTD . The contribution 
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to the energy eq.(|2T|) from the supercurrent term Scur is small for large R but becomes 
significant at denser lattices. The third term, magnetic energy Smag yields a small deviation 
of magnetic skyrmion energy from 2 at large R. 

Profile of the angle 0(p) and of the magnetic field h are depicted on Fig. 6a and 6b 
respectively. Radius of the circular cell R varies from 20A to 300A. On Fig.6a smaller value 
of R corresponds to a lower curve. Small p asymptotics of the solution up to p^ terms read: 
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where c and 6(0) are constants to be determined by numerical integration. Most of the flux 
goes through the region where the vector / is oriented upwards. In other words, the magnetic 
field is concentrated close to the center of a magnetic skyrmion. 

The value of hci{R oo) = ems{R — ^ c>o)/4 for a triplet superconductor filling the whole 
space is equal to 1/2. In physical units this result reads: 

Hci = (29) 



It is quite different from iJ^i of conventional (s-wave) superconductors where an additional 
factor logK is present. Line energy of Abrikosov vortices Ey for the present model was 



calculated numerically (beyond London approximation) in |T^. For k = 20 and 50 we 
obtain 2£y/ems ~ 3.5 and 4.4 respectively. Therefore we expect that the lower critical field 
of f/Pts is determined by magnetic skyrmions. 



C. Magnetization of the skyrmion lattice 

\ih > hci the external magnetic field enforces a definite value of magnetic flux through a 
sample. Magnetic skyrmions, being topological objects, carry quantized magnetic flux and 
their number in the sample is determined by the average magnetic induction b, similarly to 
the case of vortices. Energy of magnetic skyrmions as function of R eq.(P^) actually deter- 
mines the interaction between them. However, magnetic skyrmions, contrary to vortices, 
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are extended objects and their linear size R is also determined by the number of them in 
the sample. 

To qualitatively estimate the magnetization curve produced by "skyrmion mixed state" 
we make use of the unit cell energy obtained in the previous section. The Gibbs energy 
density of the sample of volume V = S x L, where S is the transverse area and L is 
longitudinal extension, is given in dimensionless units eq. ([lO|) by 

G{b) = -2bh=^-(2 + 5.62^) - 2bh. (30) 

The second equality follows the facts that magnetic induction b is related both to the number 
of magnetic skyrmions = Sb/2 and to the size of the magnetic skyrmion defined above 
= 4/6. Minimization of eq.(^) with respect of b yields: 

6- 0.225 ( — -ij , h>Ki = -. (31) 

Eq.(|3T|) shows that a skyrmion lattice is characterized by zero slope of magnetization 
curve at hd, in contract to the infinite slope for the magnetization curve associated with a 
vortex lattice. This circumstance provides a tool in the experimental search for the triplet 
superconductivity with approximate SO (3) symmetry. Our results agree well with the earlier 
work of Burlachkov et al who also obtained zero slop of the magnetization at h^i for a 
stripe I texture which might arise in the case of very high anisotropy of effective mass tensor 
m* (see eq.(0)). 

V. INFLUENCE OF S0(3) BREAKING TERMS 

In this section we consider influence of an 5*0(3) symmetry breaking terms on skyrmion 
lattice. List of these terms was given in section IIB eqs.(pa|-p^). The perturbations are not 
expected to affect the existence of topological solitons - just modify their energy. When 
the coefficient of a breaking term becomes of order 1, the soliton might disappear, although 
it is not necessary. We study in detail the influence of Zeeman term eq.(^). The choice 
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is motivated by our previous study of possible spontaneous vortex state in a new bulk 
perovskite superconductor Sr2Y Rui^xCu^Oq |T^. 

This compound has very unusual magnetic properties and is suspected to be a p-wave 
superconductor for the following reasons 0. At the temperature of about 60K, at which 
superconductivity sets in, these materials begin to exhibit basic ferromagnetic properties 
like hysteresis loop. Experimental observation of a positive remanence suggests existence of 
spontaneous magnetization in the absence of an external magnetic field. Exact overlap of 
superconductivity and ferromagnetism lead us to consider an isotropic triplet model eqs.(|D- 
(§) in nonunitary phase with spontaneous time reversal symmetry breaking. In this case, a 
direct spin coupling of the condensate to a magnetic field 

^S-B = ^gS-B (32) 
2m* c 

becomes relevant. In what follows this coupling will be referred to as Zeeman-like coupling 
and characterized by dimensionless parameter g. For sufficiently large values of g energetics 
of the triplet superconductor changes considerably. There exists a critical value g^ = 1 
above which the mixed state might respond on an external magnetic filed ferromagnetically 
and, on the other hand, in the presence of an external magnetic the field mixed state might 
occur even for temperatures above |jT9|. For larger Zeeman-like coupling, g > gc2 ^ ^ogK, 
vortex energy becomes negative. Spontaneous vortex phase appears at if = and exists for 
arbitrarily large magnetic field. Meissner phase, therefore, completely disappears. Vortices 
become thinner when H grows. The structure of the vortex core is markedly different from 
the usual one. 



Our analysis in |T9[ was entirely based on the simplest possible topological objects: 
vortices of the usual type. Value of k for the materials of Wu et al are estimated to be 
quite large and, consequently, vortices should be heavy compared to magnetic skyrmions. 
Spontaneously magnetized skyrmion lattice can also occur, as in the previous case of vortices 
of usual type. Values of g required to obtain spontaneous vortex state g = log k were very 
high and made the scenario questionable. This value is lowered to (7 ~ 1 for magnetic 
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skyrmion lattice. 

The free energy per unit length for a single magnetic skyrmion now has form: 



^L = j pdp 



1 fdeV sin^e /i + cose V 



/ a da\^ / a da\ ^ 
VP dp) ^ \p dp) 



The equations read: 



^„ 1^/ sin© /'2 + cos© ^ \ . f , o,\ , ^, 

e" + -9' = h 2a + sin e a' + - , (34 

P P V P / VP/ 



a" + - - ^ - a = - (1 + cose) - ^0'sin©. (35) 
P P P 2 



We use the same boundary conditions as that for the case of isolated magnetic skyrmion 
at gi = (see Sec. IV A) Calculations were performed both for positive and negative values 
of g. Plot of the energy of the magnetic skyrmion as a function of g is presented in Fig. 7. 
The characteristic feature of this dependence are a maximum near g = 0. Profiles of the 
magnetic field b{p) for different g of both signs are presented in Fig 8. Zeeman interaction 
strongly influences behavior of b{p) near the center of the magnetic skyrmion and in quite 
different manner for positive and negative g. Note, however, that as l^fl increases behavior of 
the function changes significantly in the interval of p from the origin up to only some limiting 
value, after which it remains approximately the same for different g. Thus we observe that 
nonzero g actually introduces new length scale in the problem. Changes in the profile of 
0(p) with g are less pronounced and are not displayed. 



VI. DISCUSSION 

In this paper we performed topological classification of solutions in SO (3) symmetric 
Ginzburg - Landau free energy. This model with addition of very small symmetry breaking 
terms describes heavy fermion superconductor UPts and possibly other triplet supercon- 
ductors. A new class of topological solutions in weak magnetic field carrying two units of 
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magnetic flux was identified. These solutions, magnetic skyrmions, are nonsingular (do not 
have singular core hke Abrikosov vortices). They repel each other as 1/r at distances much 
larger then magnetic penetration depth A forming relatively robust triangular lattice. At 
lattice spacings much larger then A their energy is reduced by a factor of the order of log k 
as compared to the usual Abrikosov vortex solutions and therefore dominate the magnetic 
properties for strongly type II superconductors. The lower critical magnetic field Hd = 
is reduced correspondingly by a factor 2 log k. 

Magnetization near Hd instead of sharply rising with infinite derivative increases grad- 
ually as {H — HciY- This agrees very well with the experimental results for UPts, see Fig.l 
of For fields higher then several Hd London approximation is not valid anymore since 



magnetic skyrmions will start to overlap. At distances between fluxons of order A (or at the 
field H'^i ~ Hd ■ 2 log k) one expects that ordinary Abrikosov vortices, which carry one unit 
of magnetic flux, become energetically favorable. The usual vortex picture has indeed been 
observed at high fields by Yaron et al . Curiously, our result on magnetization is similar 
to conclusions of Burlachkov et al. who investigated stripe-like (quasi one dimensional) 
spin textures in triplet superconductors. Magnetic skyrmions are quite stable objects and 
they are not destroyed by small perturbations of exact 5*0(3) symmetry of the original model 
eqs.(|]-|^). Moreover, deformed magnetic skyrmions might exist even at large deviations from 
exact 5*0(3) symmetry. We demonstrated this including Zeeman like interaction eq . (|3^) . 



Let us list below the experimental features which can allow identification of the magnetic 
skyrmions lattice. 

1. The lower critical filed is significantly smaller then usually expected. For such strongly 
type II superconductors as UPt^, Sr2RuO^ or Sr2Y Rui-xCuxOq with k ~ 50 ^ 70 the 
reduction amounts 8 times. Although Hd is expected to be very small (less then IGauss) 
it is still measurable. 

2. Magnetization above Hd-, but below crossover to Abrikosov vortex lattice H'^^ ~ 
2^ log K is markedly distinct from the usual one due to long range nature of the magnetic 
skyrmions. 
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3. Unit of flux quantization is different: 2$o- 

4. Magnetic fleld proflle is different: no exponential drop even at very sparse lattices. 

5. Superfluid density | P is almost constant throughout the mixed state. There is 
no normal cores of the fluxons. This can be tested using scanning tunneling microscopy 
technique. 

6. Due to the fact that there is no small normal core where usually dissipation and 
pinning take place, one expects that pinning effects are greatly reduced. Correspondingly 
critical current should be very small. 

7. The vortex lattice in the region around Hd can melt into so called lower fleld vortex 
liquid due to thermal fluctuations p6[. The melting of usual Abrikosov vortex lattice is easy 
even in not very strongly fluctuating superconductors because interaction between Abrikosov 
vortices is exponentially small. It is not so for magnetic skyrmions. Due to their long range 
1/r interaction the lattice is more robust and therefore no melting is expected. 
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List of Figures 

Fig. 1 Definition of angles and 0. Unit vectors /, n, m constitute a triad of perpendicular 
vectors in the spin space. ^ is the superconducting phase defined in eg. ([TB|) . 

Fig. 2 Configuration of a magnetic skyrmion with Q = —1. Solid arrows represent I field 
while "clocks" show that phase rotates twice clockwise as a round on the remote 
contour is completed. 

Fig. 3 Magnetic field of the isolated magnetic skyrmion. Distance from the center p varies 
from to a cutoff pmax with boundary conditions (9' + 6/p) \p=p^^^ = 0,b\p=p^^^ = 
imitating the infinite domain, pmax = 40A for the lowest curve and 600 for the 
uppermost one. 

Fig. 4 A fragment of the magnetic skyrmion lattice. For numerical calculations we approxi- 
mate symmetric unit cell by the disk of the same area: R = ^^^a- 

Fig. 5 Energy of the unit cell of the magnetic skyrmion lattice. Dots are numerical values for 



different R. Line is the fit of eq.(28 



Fig. 6 Numerical solution of GL equations in London approximation for a unit cell of the 
magnetic skyrmion lattice. Radius of the circular cell R varies from 20A to 300A. a) 
Angle G as a function of the distance p from the center of the cell, b) Magnetic field 
6 as a function of the distance p from the center of the cell. A smaller R corresponds 
to a lower curve. 

Fig. 7 Energy of the isolated magnetic skjrrmion as a function of dimensionless Zeeman cou- 
pling g for R/\ = 300. 

Fig. 8 Magnetic field of the isolated magnetic skyrmion as a function of distance from the 
center p for different Zeeman coupling and for the case of pmax = 300A (see Fig. 3 
caption), a) g = 0, .5, .7, .9, 1., 1.1. h) g = 0, —.5, —.7, —.9, —1., —1.1. In both cases a 
smaller \g\ corresponds to a lower at p = curve. 
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FIGURES 

Fig. 1. A. Knigavko & B. Rosenstein 
"Magnetic skyrmions and their lattices in 
triplet superconductors" 
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Fig. 2. A. Knigavko & B. Rosenstein 
"Magnetic skyrmions and their lattices in 
triplet superconductors" 
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Fig. 3. A. Knigavko & B. Rosenstein 
"Magnetic skyrmions and their lattices in 
triplet superconductors" 




Fig. 4. A. Knigavko & B. Rosenstein 
"Magnetic skyrmions and their lattices in 
triplet superconductors" 
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Fig. 5. A. Knigavko & B. Rosenstein 
"Magnetic skyrmions and their lattices in 
triplet superconductors" 




Fig. 6a. A. Knigavko & B. Rosenstein 
"Magnetic skyrmions and their lattices in 
triplet superconductors" 
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Fig. 6b. A. Knigavko & B. Rosenstein 
"Magnetic skymiions and their lattices in 
triplet superconductors" 
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Fig. 7. A. Knigavko & B. Rosenstein 
"Magnetic skyrmions and their lattices in 
triplet superconductors" 
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Fig. 8a. A. Knigavko & B. Rosenstein 
"Magnetic skyrmions and their lattices in 
triplet superconductors" 




Fig. 8b. A. Knigavko & B. Rosenstein 
"Magnetic skyrmions and their lattices in 
triplet superconductors" 
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